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Abstract 

The formation of non-relativistic collisionless shocks in laboratory with ultrahigh intensity lasers 
is studied via ab initio multi-dimensional particle-in-cell simulations. The microphysics behind 
shock formation and dissipation, and the detailed shock structure are analyzed, illustrating that 
the Weibel instability plays a crucial role in the generation of strong subequipartition magnetic 
fields that isotropize the incoming flow and lead to the formation of a collisionless shock, similarly 
to what occurs in astrophysical scenarios. The possibility of generating such collisionless shocks in 
laboratory opens the way to the direct study of the physics associated with astrophysical shocks. 
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Understanding how collisionless shocks are formed and propagate in unmagnetized plas- 
mas is of great importance to the study of many astrophysical scenarios such as gamma-ray 
bursts (GRB) afterglows, active galactic nuclei, pulsar wind nebulae, and supernova rem- 
nants l|-|3|]. The synchrotron radiation collected in astronomical observations suggests that 
these structures can generate subequipartition magnetic fields and accelerate particles to 
very high energies [4]. How the magnetic fields are generated and what their structure is, 
which dissipation mechanism is dominant, which physical processes lead to shock formation, 
and how particles are accelerated, remain open questions. Since Coulomb scattering cannot 
be responsible for the dissipation process in a collisionless shock wave, anomalous heating 
associated with particle scattering in plasma turbulence seems to be the natural explanation 

fl. 

Electromagnetic turbulence associated with the Weibel, or current filamentation, insta- 
bihty B .believed to be the .eading mechanism for shock formation in weah.y magnetized 
plasmas [7j. This instability can generate small-scale magnetic fields in counterstreaming 
plasmas which can scatter particles and isotropize the flow. Previous numerical studies of 
idealized astrophysical collisionless shock scenarios have shown, using particle-in-cell (PIC) 
codes, that Weibel instability can lead to strong filamentation, magnetic field generation, 
and shock formation js, [9] and that nonthermal particles can be accelerated in this shock 



structure 
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llj and emit synchrotron radiation 12j. The validation of these numerical 



studies for astrophysical scenarios, where in situ observations are not possible, is limited, 
since the information available from these astrophysical objects comes only from their radi- 
ation emission, which occurs at significantly larger temporal and spatial scales. 

Laboratory experiments can play a crucial role in validating theoretical and numerical 



models of astrophysical phenomena 



13l | . In the last years, a few experimental studies have 



been proposed and conducted for the generation of non-relativistic electrostatic collisionless 
shocks in laboratory with colliding laser-ablated plasmas [jjj]. However, in the case of Weibel 
mediated collisionless shocks, the conditions for shock generation in laboratory are not yet 
fully understood and no experimental evidence has been observed to date. The advent of 
high energy, high power laser systems is allowing for the exploration of extreme regimes in 
laser-plasma interactions, where strong plasma flows can be generated and the conditions 
necessary for the generation of Weibel mediated shocks in laboratory may be reached for the 
first time. Previous PIC studies of shock formation in intense laser-plasma interactions have 
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focused on the one- dimensional (ID) dynamics [15|, ll6| and, therefore, could not evaluate 
the role/impact of Weibel instability in these scenarios. 

In this Letter, we demonstrate the possiblity to generate truly Weibel mediated collision- 
less shocks in laboratory by the interaction of an ultraintense laser pulse with an overcritical 
plasma. Using ab initio multi-dimensional relativistic PIC simulations, we examine in detail 
the physics behind shock formation and propagation, from the generation of the incoming 
flow by the intense laser to the microinstabilities that lead to the generation of subequiparti- 
tion magnetic fields that isotropize the flow and lead to the formation of the shock structure. 
We show that the underlying physics is similar between these non-relativistic laser-driven 
shocks in laboratory and previously considered relativistic astrophysical shocks, illustrat- 
ing the possibility of directly studying the physical mechanisms behind these astrophysical 
scenarios in laboratory. 

In order to study the self-consistent shock formation and propagation in realistic lab- 
oratory conditions we use two-dimensional (2D) and three-dimensional (3D) simulations 
performed with OSIRIS [l7j, a fully relativistic, electromagnetic, and massively parallel PIC 
code. We simulate the interaction of an ultraintense laser pulse with a pre-ionized unmagne- 
tized electron-proton plasma. The laser is linearly polarized and has a wavelength of 1 /im. 
We have simulated different laser intensities, ranging from 10 20 — 10 22 Wcm -2 , correspond- 
ing to a normalized laser vector potential a$ ~ 10 — 100, and different plasma densities, 
ranging from 10 — 100 n c , where n c is the critical density for 1 /zm light. The initial plasma 
temperature is 1 keV. The computational domain is typically 80 c/u P i in the longitudinal 
direction and 18 c/u P i in the transverse direction, with c/u P i = c{A^Z 2 e 2 n p /rrii)~ 1 ^ 2 the ion 
skin depth for a plasma density n p and ion mass m,; e is the elementary charge, Z is the 
charge state, and c the speed of light in vacuum. The system is numerically resolved with 2 
cells per c/u pe in both directions and uses 64 particles per cell for each species, for a typical 
total of 10 9 particles. We use cubic particle shapes and current smoothing with compensa- 
tion. Larger transverse box sizes, higher resolution, and higher number of particles per cell, 
were tested, showing overall result convergence. We note that these are the first full-scale 
simulations of unmagnetized electron-ion shocks (previous simulations used either positrons 
or ions with reduced mass ratios). 

As the intense laser hits the overdense target it acts like a piston, pushing the front of the 
target as a massive and uniform flow but also generating a population of fast /hot electrons. 
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Fast electrons, which typically have a density on the order of the critical density, n c , and 
a relativistic factor 70 ~ \Jl + of, are not affected by the proton response and propagate 
through the target. A cold return current is set up in order to balance the incoming fast 
electron flux. The two counterstreaming flows go Weibel unstable, similarly to what is 
believed to occur in astrophysics 0-Q- Even in the case where the laser generated incoming 
flow is too hot to filament by itself in the background plasma, the return current is cold 
and therefore filaments. Following the usual procedure for the calculation of the dispersion 
relation for purely transverse modes 18|, and in the limit of a hot and rarefied electron flow 
counterstreaming with a cold and dense slowly drifting electron background, the maximum 
electron Weibel instability growth rate can be shown to be simply Twe — ifirl \Jlr)u pe - l 
where /3 r is the normalized velocity of the returning electrons and 7 r their Lorentz factor. 
In order to establish current neutrality (3 r ~ n c /n p , yielding T We = 0.01 — 0.1 u pe for 
the parameters in study, and leading to the generation of strong magnetic fields in a few 
100 u) pi }. The electron instability saturates when the magnetic energy density in the Weibel 
filaments becomes comparable to the energy density in the fast electron flow, leading to a 
saturation magnetic field amplitude B sat ~ (87raon c m e c 2 ) 1//2 , for ultrahigh laser intensities 
(a ^> 1). The magnetic fields associated with Weibel/current filamentation instability of 
the fast electron flow isotropize the incoming non-relativistic electron-proton flow, leading 
to a strong compression and to the formation of a shock, defined as the density compression 
that propagates away from the laser-plasma interface. The shock speed is determined by the 
slowdown and mass/pressure build-up associated with this more massive flow, and therefore 
it is non-relativistic. Once the shock is formed, the particles that escape the shock from the 
downstream still provide the generation of the counterstreaming cold return current in the 
upstream, which allows for continuous filamentation in the upstream region. 

For the collisionless shock to be formed in the unmagnetized plasma, it is required that 
the piston (downstream) velocity exceeds the ion sound speed, eg = (ZkbTe/mi) 1 ^ , where 
kf, is the Boltzmann constant, and T e is the bulk electron temperature, which is typically a 
fraction a ~ 1/3 — 1/2 of the ponderomotive temperature before the shock is formed 19], 
kbT e ~ aaom e c 2 . The downstream velocity, v^, can be estimated by equating the momentum 
flux of the incoming mass flow with the laser light pressure, yielding a normalized velocity 



of /3d = Vd/c = a/ (n c /2n p )(Zm e /mi)al [20j. The condition for shock formation, Vd > C5, is 
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then given as a function of laser and plasma parameters by 



n.,. 



a > 2a^. (1) 

n c 

Fig. [1] illustrates the main features of shock formation for a typical simulation where we 
have used a laser intensity of 5 x 10 21 Wcm -2 (a = 60) and a plasma density of 50 n c . A 
strong compression is observed in the downstream (Figs. [T] a and[T]b), between the laser- 
plasma interface and the shock front, and strong filamentation in the upstream region. The 
magnetic field illustrates similar filamentary structures (Fig. d]c). The filaments of density 
and magnetic field are not stationary in front of the shock. The filaments size evolves 
from the electron skin depth c/u pe , far upstream, to the ion skin depth c/u P i, close to the 
shock, and are then frozen behind the shock front. The shock transition is about 1-2 c/u P i 
thick, at early times, which is of the order of the ion Larmor radius, and corresponds to 
a peak in the magnetic energy (Fig. [Ud). At later times, the thickness of the magnetic 
field peak continuously increases towards the downstream region, as observed in previous 



astrophysical configurations 21], reaching values of the order of 10 c/u P i for our largest 
interaction times, 500 c^C . The transversely averaged magnetic field energy density reaches 
12 % of equipartition with the upstream kinetic energy density (measured in the downstream 
rest frame), i.e. e B = (B 2 /8n) / 1 (n p niiC 2 (^d — 1)) — 0.12 (Fig. Q]d). Note that for non- 
relativistic shock velocities, the upstream velocity in the downstream frame is simply the 
opposite of the downstream velocity in the laboratory (upstream) frame, which in this case 
is measured to be Vd ~ 0.13 c. Locally, the magnetic field energy density at the shock front 
can reach equipartition with the upstream, e B _ * 1. These values are fn.ly consistent with 
previous simulations of Weibel mediated relativistic shocks in astrophysical scenarios 8-11 . 
indicating that the underlying physical mechanisms are similar, even if in the present scenario 
the shock velocity is non-relativistic. An important difference is the well defined longitudinal 
electric field observed at the shock front in our simulations (Fig. [T]e). This is associated 
with the fact that downstream electrons are significantly hotter than ions, since the laser 
predominantly heats electrons. This was not observed in previous simulations of relativistic 
counterstreaming plasmas, as both electron and ion flows are initialized completely cold, 
which will hardly be the case in a laser-driven laboratory configuration. Although the energy 
associated with this electric field is relatively small, e# = (E 2 /8ir)/(n p miC 2 ('jd — 1)) — 0.025 
(Fig. H]f), the field is able to reflect a fraction of the upstream ion population (10-15%). 



As these reflected ions counterstream with the background plasma, they will lead to an 
enhancement of the magnetic fields in the ion time scales due to Weibel instability in the 
ions. In the limit of cold ions, the maximum ion Weibel growth rate is given by T^i — 
Pby/ ( n b/ n P ) I (lb m i) u pe ~ 2/3 s h a/ ( n b /n p )/mi u pe , where n b is the density of the reflected 
ion beam, which moves with twice the shock velocity. We can indeed see that in the foot 
of the shock (region where the reflected ions are present) the Weibel magnetic fields are 
stronger than in the remaining upstream region (Fig. [T]c). At this point the instability 
becomes similar to more conventional scenarios with two counterstreaming plasma flows, 
but where the electrons are relativistic and ions are non-relativistic. It should be noted that 
an electrostatic ion-ion instability associated with the reflected ions has been previously 
identified in electrostatic shocks, strongly affecting its structure 22]. It can be shown, 
following the usual procedure for the calculation of the dispersion relation for electrostatic 



modes 



221 ]. that in the limit of cold ion flows, the maximum growth rate of the ion- ion 



electrostatic instability is Tei — \J {nb / n p ) / (8mj) u pe . This instability tends to dominate 
over the Weibel in the case of low shock velocities ((3 s h < 0.1); however, for the large shock 
velocities reached in our proposed setup ((3 s h > 0.1) the Weibel instability dominates, further 
amplifying the magnetic fields and isotropizing the incoming flow. Even in the case where 
the initial maximum growth rates are comparable, the Weibel instability tends to dominate 
as the ion heating caused by the Weibel instability quickly shuts down the electrostatic 
instability. We have simulated the propagation of an ion flow in a plasma background to 
confirm this (not shown here). For the densities and velocities associated with our setup, 
and for different ion temperatures (1 eV - 1 keV), we consistently observe that the Weibel 
instability dominates over the electrostatic ion-ion instability, whereas in the case of lower 
ion velocities (/3b < 0.2) the electrostatic instability dominates. The condition [3 s h > 0.1 
effectively defines a lower limit for the laser intensity required to drive Weibel mediated 
shocks in this configuration. The detailed comparison of these instabilities in scenarios with 
high velocity ion flows will be given elsewhere. 

The particle spectrum at different longitudinal positions is highly modified by the shock 
structure. In the downstream region (Fig. Eg), we observe a two-temperature electron spec- 
trum from the laser acceleration, which can be reasonably fitted to a sum of 2D Maxwellians 
of the form f(j) = a!7exp(— 7/A71) + a 2 7exp(— 7/A72), with normalizations a\ and a 2 , 
A71 = 13, and A7 2 = 58, which is close to the expected laser induced ponderomotive tem- 
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perature of A7 = a/1 + ag ~ 60 20]. The bulk electron temperature, A71, which initially is 
a fraction a of the ponderomotive temperature, changes as the shock is formed and most of 
the particles are trapped behind it, leading to an equipartition between the electron thermal 
energy and the ion fluid energy, i.e. A71 = (7^ — l)mj/m e . The most energetic electron 
population crosses the shock front into the upstream and remains relatively unchanged as 
evidenced by the spectra at the foot of the shock (Fig. [T]h) and in front of it (Fig. [TJi). The 
ions are heated up in the downstream (Fig. (Hg), whereas in the upstream we observe a cold 
ion Maxwellian background and the presence, in the foot region, of a small ion population 
that has been reflected at the shock front (Fig. [T]h). In order to compare directly our 
results to previous astrophysical simulations, we have transformed the different quantities 
into the downstream frame, in order to calculate the energy balance between electrons and 
ions at the different regions. As the filaments merge in the vicinity of the shock front and 
the magnetic fields of the filaments coalesce, the energy in the fields is converted back to 
thermal energy of the ions in the downstream region, thus leading to an effective heating of 
the ions and electrons. We observe that in the upstream region, the initial ion flow moving 
towards the shock loses 20-25% of its energy for electron heating, and another 20-25% goes 
into ion heating (thermal energy) in the downstream. These numbers are comparable, but 
lower than the observed values in idealized astrophysical configurations {23)]. 

The observed features illustrate the collisionless nature of this shock structure, which has 
a transition region that is significantly smaller than the typical particle-particle scattering 
mean free path. The typical temperature and density behind the shock front is observed to 
be Ti ~ 10—100 keV, T e ~ 1 — 10 MeV and n p ~ 30— 300n c , corresponding to a typical proton 
(electron) Coulomb mean free path on the order of 10 5 (10 6 ) c/u P i. In front of the shock the 
temperature can be significantly lower. In our simulations, with keV pre-ionized plasmas 
with densities of the order of a few tens of n c the proton (electron) Coulomb mean free path 
is still ~ hundreds (few) c/u P i, and therefore no significant collisional effects are expected. 
We redid our simulations with a relativistic Monte Carlo Coulomb collisional operator 24] 
in order to accurately model the effect of collisions and no significant differences have been 
observed. In this collisionless shock structure dissipation is provided by particle scattering 
in the self-generated Weibel turbulence. In other setups, where cold (10 - 100 eV), solid 
density (> 100n c ) plasmas might be produced, v ei jijj pe ~ 0.1 — 1 and therefore collisions 
between fast particles and the cold background plasma should be taken into account. 



7 



The generated shock structure has a well defined velocity and density jump. The shock 



hydro dynamic jump conditions 



25] predict a density jump with a downstream to upstream 



density ratio n 2 /n 1 = (T ad 'y d + l)/(T ad — 1) and normalized velocity of f3 s h = v s h/c = 
[(1 + T ad -f d )(-fj - l) 1/2 )]/[l +J d + Tad(lj ~ 1)], where T ad is the adiabatic index and j d is the 
Lorentz factor of the downstream in the frame of the upstream (which for this case is the 
laboratory frame, since the upstream is approximately at rest). In the non-relativistic limit, 
where (3 d <C 1 and 7^ ~ 1, the density jump can be written as n 2 /n 1 = (T ad + l)/(T ad — 1) 
and the shock velocity as 

a I n c Zm e 

provided that the plasma is always opaque to the incoming light, i.e. that the downstream 
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density is higher than the relativistic critical density n p (T ad + l)/(T ad — 1) > n c y\ + a . 
More general jump conditions, where the adiabatic index depends on the details of the 
distribution function of the downstream, can be used to infer the exact jump conditions 



corrected for distributions functions that deviate from Maxwel 
the different plasma species and will be presented elsewhere 



ian and/or are different for 



26j | . For the case of Fig. [fl we 



observe a density jump of ~ 3.1 and a normalized shock velocity of ~ 0.19 (Fig. [2]), which 
are in very good agreement with the derived theoretical values of 3 and 0.2 respectively, for 
an adiabatic index of 2, appropriate for a 2D gas. 

In order to understand if our 2D simulations can capture the relevant 3D physics, we have 
performed 3D simulations for the same parameters of Fig. [Q A similar shock structure has 
been obtained with a normalized shock velocity of 0.18, which is consistent with Eq. (j2J) for 
an adiabatic index of 5/3, appropriate for a 3D gas. Conversely, ID simulations of the same 
scenario cannot capture the relevant physics, in particular the current filamentation/Weibel 
instability, and therefore yielded different results. 

As evidenced by the shock jump conditions, the properties of Weibel mediated collisionless 
shocks generated in laboratory in the configuration here proposed can be controlled by 
tuning the laser and plasma parameters. In particular the shock velocity can be controlled 
by adjusting the laser intensity, /, the plasma density, n p , and/or the target composition (ion 
mass, m^, and charge state, Z) according to Eq. (j2J), affecting proportionally the energy of 
the ions reflected by the shock, e ion [MeV] ~ 74.2Z/[10 21 Wcm- 2 ]/n p [10 22 cm 1 , which can be 
measured experimentally. This tunability is illustrated in Fig. [31 where we can see that the 
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shock velocity obtained for different simulations with different laser and plasma parameters 
is in very good agreement with Eq. ([2]). Furthermore, the laser polarization can also be 
varied (from linear to circular) in order to control the laser electron heating and to study 
its influence in the shock properties. 

The laser parameters required to study the generation of Weibel mediated collisionless 
shocks in laboratory can be estimated based on our results. The laser duration, To, should 
be significantly longer than the shock formation time, tf, which is of the order of the time 
it takes the piston to push the plasma a shock thickness forward, tf ~ (/^Wpi) -1 , i.e. 
r > lO^dUpiy 1 or r [ps] > 0.5^- (/ [10 21 Wcm _2 ]Z 3 ) , where m p is the proton mass. 
The laser spot size should be much larger than the shock thickness in order to have a 
stable shock front, W > Wc/ou pi or W [/H > 10 {^) ^ Combining 

the two previous conditions, the laser energy required to provide a stable shock structure 
is e [kJ] ~ 1.5 ^o[io 21 Wcm-^ y/ 2 ^x [fim}) 2 (55^) ~\ We note that high-Z materials 

allow the use of lower energy lasers in order to drive a stable shock structure. Current 
and near-future laser systems can deliver the required picosecond 100 J - kJ pulses, either 
directly or using a recently proposed Raman amplification scheme to convert nanosecond 
kJ pulses into picosecond kJ pulses 



281 ] . We also note that the simulations here shown 
are valid for different laser wavelengths and plasma densities, provided that the laser to 
plasma frequency ratio is kept constant and that collisional effects can still be neglected. 
For instance, C02 lasers interacting with gas jets, which can easily provide densities above 



the critical density 



shock structures [291 ] . 



or 10 /iin laser wavelength, could be used in order to generate similar 



The versatility of the setup here presented should also allow for the generation of lab- 
oratory shocks in different astrophysical relevant regimes. Namely, the use of non-uniform 
targets would allow for the study of shock formation and propagation in a clumpy medium 
where magnetic energy production by macroscopic turbulence may become important 30] , 
and the use of strong, externally induced, uniform magnetic fields 31[ can be envisioned to 
study shock formation and propagation in plasmas with variable magnetization. 

In conclusion, we have studied the generation of Weibel mediated collisionless shocks 
in unmagnetized plasmas, driven by the interaction of an ultraintense laser pulse with an 
overcritical plasma. We have shown that in realistic laboratory conditions the plasma flow 
generated by this interaction can lead to the generation of subequipartition magnetic fields, 
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due to Weibel/current filamentation instability, that isotropize the flow and generate a 
shock. The shock structure and its properties, here shown for the first time for realistic ion 
to electron mass ratios, are similar to previously simulated low mass ratio Weibel mediated 
collisionless shocks in idealized astrophysical scenarios. Our results illustrate the possibility 
of studying for the first time in laboratory the physics behind the formation and propagation 
of Weibel mediated collisionless shocks in unmagnetized plasmas, which would allow for a 
better understanding of the role of these structures in nonthermal particle acceleration and 
emission of synchrotron radiation in astrophysical scenarios. 
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FIG. 1: Steady state structure of a collisionless shock generated through the irradiation of an 



overcritical unmagnetized plasma by an ultraintense laser pulse after 385 fs (113 u> ^) of interaction. 



pi 



a) Density structure normalized to the unperturbed upstream density, b) Transversely averaged 
plasma density, c) Magnetic field in the direction outside the simulation plane, d) Transversely 
averaged magnetic energy, e) Longitudinal electric field, f) Transversely averaged electric energy, g- 
i) Electron and ion (scaled up by the mass ratio mi/m e ) spectrum at three different slices (positions 
marked by arrows). Red: electrons; blue: ions; green: fit to a sum of two 2D Maxwellians. The 
transversely averaged quantities in b), d), and f) are also shown for interaction times of 805 fs (250 
o;^ 1 ) (dashed) and 1134 fs (353 w^ 1 ) (dotted). 
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t [1 / co pi ] 

FIG. 2: Temporal evolution of the transversely averaged density profile illustrating the normalized 
velocity of the downstream /J^, the shock (3^, and the ions that are reflected at the shock front j3 r . 




n p /n c 

FIG. 3: Shock velocity as a function of the plasma density and the laser normalized vector potential. 
Numerical values on the plot indicate the shock velocity measured in simulations. The parameter 
range has an upper limit defined by the condition for relativistic opacity of the downstream plasma 
(solid) and lower limits defined by the condition for the Weibel instability to dominate over the 
electrostatic instability, ji^ > 0.1 (dotted), and by Eq. (P), with a = 1/3 (dashed). 
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